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Hydrogen-helium mixtures at conditions of Jupiter's interior are studied with first-principles com- 
puter simulations. The resulting equation of state (EOS) implies that Jupiter possesses a central 
core of 14 - 18 Earth masses of heavier elements, a result that supports core accretion as standard 
model for the formation of hydrogen-rich giant planets. Our nominal model has about 2 Earth 
masses of planetary ices in the H-He-rich mantle, a result that is, within modeling errors, consis- 
tent with abundances measured by the 1995 Galileo Entry Probe mission (equivalent to about 5 
Earth masses of planetary ices when extrapolated to the mantle), suggesting that the composition 
found by the probe may be representative of the entire planet. Interior models derived from this 
first-principles EOS do not give a match to Jupiter's gravity moment J4 unless one invokes interior 
differential rotation, implying that jovian interior dynamics has an observable effect on the measured 
gravity field. 



I. INTRODUCTION 



The discovery of over two hundred extrasolar planets 
resembling the giant planet Jupiter in mass^ and compo- 
sition^ has raised fundamental questions about the origin 
and inner structure of these bodies. Establishing the exis- 
tence of a dense core in Jupiter is vital for understanding 
the key processes of planetary formation. An observed 
correlation between the metallicity of the parent star and 
the likelihood of giant planets orbiting it^*^ may indicate 
that an initial core aggregated from solid planetesimals 
triggers the gravitational collapse of nebular gases to 
form giant planets^i, or alternatively, increased metal- 
licity enhances the direct collapse of giant planets from 
the nebula--. Current planetary models^i^^ predict only 
a very small core between and 7 Earth masses (Me) for 
Jupiter, lending support to core-accretion theories with 
comparatively small cores2^, late-stage core erosion sce- 
nariosi^, or suggesting that jovian planets are able to 
form directly from gases without a triggering core^. Us- 
ing first-principles simulations for hydrogen-helium mix- 
tures at high pressure, we show here that Jupiter pos- 
sesses a significant central core of 16±2 Me of heavier el- 
ements, a result that supports core accretion as the stan- 
dard model for formation of hydrogen-rich giant planets. 
In our model, Jupiter's mantle, defined as the H-He-rich 
outer layers, shows no enrichment of heavier elements, in- 
dicating that compositional measurements from the 1995 
Galileo Entry Probe mission may be representative of 
most of the planet. The derived interior models also pro- 



vide the first evidence that Jupiter's interior does not ro- 
tate as a solid body because no match to Jupiter's gravity 
moment J4 is obtained unless one invokes interior differ- 
ential rotation, implying that jovian interior dynamics 
have an observable effect on the measured gravity field. 

While laboratory techniques cannot yet probe deep 
into Jupiter's interior (reaching pressures P ^ 100-1000 
GPa and temperatures T ~ 10000 K), advances in first- 
principles computer simulation techniques have made it 
possible to characterize the equation of state (EOS) of 
hydrogen-helium mixtures for the entire planet. The de- 
termination of the EOS, in combination with the mod- 
eling reported here, enables us to more precisely specify 
Jupiter's mantle metallicity and core mass. 



II. SIMULATION METHODS 

All density functional molecular dynamics (DFT-MD) 
simulations were performed with either the CPMD 
code (") using local TrouUier-Martins norm-conserving 
pseudopotentials or with the Vienna ah initio simulation 
packageiS. using the projector augmented-wave method^. 
The nuclei were propagated using Born-Oppenheimer 
molecular dynamics with forces derived from the elec- 
tronic ground state. Simulations with thermally excited 
electronic states did not yield any significant correction 
to the EOS in Jupiter's interior. Exchange-correlation ef- 
fects were described by the generalized gradient approx- 
imation^. The electronic wavefunctions were expanded 
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FIG. 1: Isochores derived from DFT-MD simulations of H-He 
mixtures (circles, Y=0.2466) and pure hydrogen (triangles). 
Results for mixtures predict a positive Griineisen parameter, 
[dP /dT)\v > 0, along Jupiter's isentrope (dashed line) pre- 
dicting the planet to be isentropic and fully convective. 



in a plane-wave basis with a 35-50 Hartree energy cut- 
off. The simulations were run for 2 picoseconds with 
time steps ranging between 0.2 and 0.8 femtoseconds. We 
performed well over a hundred separate DFT-MD simu- 
lations on a non-uniform grid in density and temperature 
ranging from /9=220 - 6006 kgm^^ j^^^j 599 - 2OOOO K. 
At lower densities, we used classical Monte Carlo simu- 
lations with fitted potentials^i, which are in good very 
agreement with the SC EOS. 

All simulations were performed with 110 hydrogen and 
9 helium atoms in periodic boundary conditions. The ini- 
tial simulations were performed with F point sampling of 
the Brillouin zone only. Further simulations with uni- 
form grids of up to 4x4x4 k-points exhibited a correction 
to the pressure in the metallic regime. This correction 
slowly increases with density, eventually reaching —1.6% 
at conditions of Jupiter's core. Finite size effects were 
further analyzed by Vorberger et al.^. The isentropes 
were derived from a thermodynamically consistent fit to 
the free energy using an approach that was described in 
detail in Refi^. 

DFT-MD predicts a continuous transition from the 
molecular to a dissociated regime in fluid hydrogen as 
function of pressur o^^i^^ , which is consistent with quan- 
tum Monte Carlo simulations^. The dissociation transi- 
tion in dense hydrogen is driven by the increasing over- 
lap of molecular orbitals. At sufficiently high compres- 
sion, Pauli exclusion effects trigger a delocalization of the 
electronic charge--, the band gap close o^^'^^ , the conduc- 
tivity increases, the protons are no longer paired, and 
eventually the system assumes a metallic stated. As 
the temperature of the fluid increases, this transition oc- 
curs at lower pressures due to stronger collisions. When 
helium is added to dense hydrogen, it localizes the elec- 
tronic charge because of its stronger binding. It dilutes 
the hydrogen subsystem by reducing overlap of molecu- 
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FIG. 2: Comparison of the DFT-MD EOS with the SC and 
Sesame models. Three isotherms for pure hydrogen are shown 
in the metallic regime at high pressure. P is scaled by the 
volume per electron to the power | to remove most of the 
density dependence. 

lar orbitals, increasing the stability of hydrogen molecules 
compared to pure hydrogen at the same P and T~. 

While the DFT-MD EOS has a positive compressibility 
for all P and T, such that {dP/dp)T > (where p is the 
mass density) , the method predicts a negative Griineisen 
parameter, equivalent to {dP/dT)p < 0, for the dissocia- 
tion region of pure hydrogen- A negative Griineisen 
parameter would introduce a convection barrier into 
Jupiter's mantle. However, Figure [T] demonstrates that 
in a H-He mixture, the region where {dP/dT)p is neg- 
ative occurs at higher pressures and significantly lower 
temperatures than are present in Jupiter's interior. Con- 
sequently, Jupiter's mantle is predicted to be fully con- 
vective and isentropic. 

III. INTERIOR MODELS 

Were Jupiter of exactly solar composition and all its 
magnesium-silicates and iron together with the abundant 
hydrides CH4, NH3, and H2O concentrated in a dense 
core, considering recent reductions in the solar C and O 
abundances, its core mass would only comprise about 3 
M b^^i^^'^° . The core mass thought to trigger nebular 
collapse is several times larger^. Detection of a central 
core with a resolution of a few Me, only 1% of Jupiter's 
total mass of 318 M^, places stringent demands on the 
accuracy of the hydrogen-helium EOS. We show in this 
Letter that an EOS based on density functional molecular 
dynamics (DFT-MD) differs from the previous Saumon- 
Chabrier-Van Horn (SC) EOS^''^^ by considerably more 
than one percent. 

Figure [2] compares the DFT-MD results for pure H 
with the SC EOS. The SC model relies on analytical tech- 
niques that describe hydrogen as an ensemble of stable 
molecules, atoms, free electrons, and protons. Approxi- 
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FIG. 3: Comparison of the resulting Jupiter models. Solid 
curves show mass density (left-hand ordinate) and dashed 
curves show temperature T (right-hand ordinate), as a func- 
tion of pressure P. In contrast to Fig. [21 the DFT-MD (red 
hne) and the SC models have differing chemical compositions 
and temperatures at each pressure. The outermost layers at 
P < 100 GPa make about 90 percent of the contribution to 
the gravity moment J4. The black curves show two models 
based on the smoothed and discontinuous versions of the SC 
EOS. 



mations are made to characterize their interactions. Con- 
versely, with DFT-MD one simulates a fully interacting 
quantum system of over a hundred electrons and nuclei. 
The isotherms shown in Fig. [5] specify the EOS in the in- 
ner regions of Jupiter where hydrogen is metallic. Due to 
the lack of experimental data at this density, predictions 
from chemical models including the Sesame database^S 
vary substantially. The DFT-MD method obtains higher 
densities than the SC EOS for much of the pressure range 
characteristic of the interior of Jupiter. This has great 
influence on the predicted core mass. Moreover, DFT- 
MD predicts a continuous molecular-to-metallic transi- 
tion while the original SC EOS model predicts this tran- 
sition to be of first order, which introduced the possibil- 
ity of having different chemical compositions in the inner 
and outer layers of Jupiter. Guillot et al.— showed that 
having flexibility to redistribute materials allows one to 
make a coreless model for Jupiter while typical models 
with constant chemical composition in the mantle predict 
a core of about 7 Me when the SC EOS is used. Con- 
versely, our DFT-MD EOS predicts Jupiter's mantle to 
be isentropic, fully convective, and of constant chemical 
composition. As we show below, the mantle metallicity 
is so low that compositional gradients would be unlikely 
to change this conclusion. 

Jupiter isentropes from different models are compared 
in Fig. [3] The chemical composition corresponding to 
each isentrope was determined by adding ~ 0.7 mass per- 
cent of CH4, NH3, and H2O (collectively called "plane- 
tary ices">2S,), using ideal mixing, to the mixture of H and 
He considered here. The adequacy of the ideal-mixing ap- 



proximation was verified by several H-He simulations at 
relevant pressures and temperatures with "guest atoms" 
of C, N or O. The DFT-MD EOS was derived with a 
single H-He composition, corresponding to 9 He atoms 
per 110 H atoms, or a He mass fraction Y ~ 0.2466, to 
be contrasted with the Jovian atmospheric He mass frac- 
tion measured by the Galileo entry probe, Y = 0.2315 
± 0.0061^. The DFT-MD isentrope was then perturbed 
to the probe He abundance by using ideal mixing with a 
pure He density at the same T and P, synthesized from 
separate DFT-MD simulations for pure helium^'^. As can 
be seen in Fig. O DFT-MD implies lower temperatures 
in the deep jovian interior in isentropic models, further 
contributing to generally increased density with respect 
to SC models. The high density requires us to adopt Y 
= 0.2315 throughout the mantle in order to allow a non- 
negative planetary ice fraction. Figure [3] compares the 
original discontinuous SC EOS as well as the smoothed 
SC EOS version that we use in Figs.Oand discuss below. 

The DFT-MD EOS yields a significant revision of the 
interior structure of Jupiter. The resulting model pre- 
dicts a large central core of 16 Earth masses and a low 
metallicity for the mantle while models based on the SC 
EOS imply a small core and a much higher metallicity 
in the mantle, comprising some tens of Earth masses of 
elements heavier than H and Hei^i^^. DFT-MD based 
models predict that the planetary ices were primarily in- 
corporated into a massive solid core and depleted gaseous 
nebula. Conversely, models based on the SC EOS imply 
that planetary ices were largely accumulated along with 
the nebular hydrogen and helium when Jupiter formed. 
Figure [3] shows that the Jupiter model based on the DFT- 
MD EOS is about 5% denser than the model based on 
the smoothed SC EOS over the pressure range 30 GPa 
< P < 300 GPa, which spans about 30% of the mass of 
Jupiter, while it is a few percent less dense at deeper lay- 
ers. The Gahleo Probe found abundances of CH4, NH3, 
and H2O (at the deepest level) corresponding to at most 
about 5 Me of planetary ices if extrapolated to the entire 
mantle. Our nominal model has about 2 Me of planetary 
ices in the mantle. Given the cumulative uncertainties of 
the DFT-MD EOS and the planetary ices EOS, the nom- 
inal model is consistent with such an extrapolation of the 
Galileo measurements. 

If Jupiter formed in a region of the solar nebula where 
temperatures were too low for H2O to be in the gas phase, 
it would follow that a primordial core would include most 
of the H2O in solid form, and that accreted nebular gas 
would accordingly be depleted in H2O. This is consistent 
with the last model in Table[Tl We estimate that the non- 
H-He component of Jupiter's mantle comprises 2 ± 2 Mb, 
and Jupiter's dense core comprises 16 ± 2 Me. The er- 
ror bars include systematic and statistical uncertainties 
in the DFT-MD EOS. The dense core is modeled as rock 
and H2O in solar proportions, but could for example in- 
clude ~ 5Me of sedimented He for an initial Y = 0.25. 
Simulation parameters were chosen in order to reach an 
accuracy of 1% within the DFT method. A uniform 1% 
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FIG. 4: Jupiter's zonal wind speeds as a function of the 
distance from the rotation axis. The oscillating curve shows 
an average of northern and southern hemisphere surface zonal 
winds as determined by cloud motions; the heavy solid curve 
shows a fit of the surface winds with an eighth-order poly- 
nomial in I. The dashed curve shows our preferred rotation 
model that provides a match to Jupiter's low-degree gravita- 
tional moments with the DFT-MD EOS. The upper abscissa 
gives the relative mass enclosed within a cylinder of radius I. 



our model, we cannot match J4 by distributing helium 
unevenly in an upper and lower mantle layer, which is 
one major difference compared to earlier models e.gi^. 
Because the mantle metallicity is so low, neither can we 
redistribute mass by invoking composition gradients. It 
is more plausible that Jupiter's J4 is affected to a mea- 
surable extent by zonal winds in this pressure range. We 
bring the calculated J4 into agreement without apprecia- 
bly changing any other parameters of the model by as- 
suming the zonal wind profile shown in Fig. [¥] (preferred 
model in Table |T| . The shallower part of the wind pro- 
file could be adjusted further to force better agreement 
with Jq, but such a refinement lies outside the scope of 
this paper. In that regard, we note a recent paper— that 
models deep zonal flows in Saturn by fitting gravity data. 

To ensure that numerical errors do not play any role in 
the mismatch of J4, we checked our calculations with two 
independent theories for rotationally-distorted models of 
Jupiter. The standard theoryS^ solves iteratively for the 
shape of level surfaces to third order in the rotational per- 
turbation, sufficient to calculate J2 to a relative precision 
of 10"^, J4 to about 10^2, and Jq to about 10"^ The 
results agree well with the more accurate SCF method. 



alteration of the pressure changes the core mass by 1 A/e- 
An earlier Jupiter model based on a completely different 
set of DFT-MD EOS results for pure hydrogei*^ yielded 
a core mass that differed by 2 Me- We use this difference 
to estimate the uncertainties of our predictions. 

An acceptable model of Jupiter's interior must match 
the observed constraints of the planet's mass, equato- 
rial radius at a standard pressure of 1 bar, and the ob- 
served multipole moments J2, J4, Jq of the gravity field 
normalized to the equatorial radius. Jupiter's multipole 
moments are primarily a response to the planet's ro- 
tation, and the higher-order moments can be strongly 
affected by nonuniform rotation. Traditional models of 
Jupiter's interior structure and external gravity have as- 
sumed solid-body rotation equal to the rotation rate of 
the magnetic field (the corresponding rotation period is 
9:55:29.7 hours, stable over decadesSA). Our model calcu- 
lations find the axially-symmetric mass density through 
a self-consistent field (SCF), calculated to fourth order 
in the rotational distortioni^ii^. We incorporate differ- 
ential rotation on cylinders to the same order as the SCF. 
Table U shows that models with solid-body rotation and 
the DFT-MD EOS predict a | J4I that is seven standard 
deviations larger than the observed | J4 1 . The error bar of 
J4 has decreased by about a factor three^^ with respect 
to earlier determinations^. As predicted by theory^, re- 
duced IJ4I is associated with a more rapid decrease of 
density with radius in the pressure range 1 to 100 GPa. 
However, we have no physical justification for making 
ad hoc modifications to the P(p) relation in this pres- 
sure range. Since we have a fully convective mantle in 



It has been showni^ that high-order Jupiter gravity 
harmonics are sensitive to interior dynamics in the out- 
ermost planetary layers at pressures less than about 10 
GPa. There is no significant uncertainty in the hydrogen 
EOS at these pressures, which are accessible to experi- 
ment^^. This is why we investigated whether the DFT- 
MD EOS, coupled with plausible models of Jovian man- 
tic dynamics, could reproduce the observed Jovian grav- 
ity field. According to the Poincare-Wavre theorem^, 
the requirement that there be a single pressure-density 
relation (barotrope) within Jupiter is a necessary and 
sufficient condition that deep windspeeds be constant 
on cylindrical surfaces of constant /. Mapping the ob- 
served windspeeds onto such cylinders and then fitting 
the implied centrifugal potential with a polynomial to 
degree eight in ^, we obtain the heavy curve in Fig. U) 
When this fit is incorporated in the SCF calculations, 
we find insignificant shifts in the predicted values of J„ 
(n = 2, 4, 6), although the winds could affect higher-order 
Jupiter gravity components'^. The assumed retrograde 
deep zonal wind has negligible effect on Jupiter's to- 
tal spin angular momentum, reducing the latter by only 
0.06% with respect to a windless model. 

New observational data concerning these theoretical 
predictions are expected from the NASA mission Juno. 
Juno is a low-periapse Jupiter orbiter that will return 
unprecedented data on Jupiter's high-order gravitational 
and magnetic fields during 2016. It may present the first 
direct evidence of deep interior zonal flows in Jupiter pro- 
posed here. 
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TABLE I: Constraints on Jupiter interior structure and par- 
rameters of models based on DFT-MD EOS. Predicted con- 
ditions at the core-mantle boundary (CMB) are listed last. 



Model 


Equatorial 

radius 

(km) 


J2X10'' J4X10*' JsXlO" 

Gravity moments 


Core 
mass 
(Me) 


Mantle 

ice 

(Me) 


r(K) 

at 

CMB 


P (GPa) 
at 

CMB 


Observed 


71492 


14696.43 
± 0.21 


-587.14 
± 1.68 


34.25 
± 5.22 










Solid-body rotation 


71493 


14718 


-620 


37.5 


16.7 


2.1 


12500 


3800 


Average surface winds 


71493 


14737 


-623 


37.3 


16.7 


2.1 


12500 


3870 


Preferred model: deep winds 


71492 


14697.14 


-586.3 


23.9 


16.2 


2.1 


12500 


3800 



